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Abstract 

Galaxies are lighthouses that sit atop peaks in the density field. There is good observational evidence that 
these lighthouses do not provide a uniform description of the distribution of dark matter. 
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1 Let there be Light 

In the beginning there was dark matter and gas but there 
were no stars. Today some of the gas that was once dis- 
persed has transformed into stars, and we have light. In 
places the stars are young, hot, and bright while elsewhere 
the stars are old, cool and faint. In places the stars have 
been scattered by violent collisions. In places the gas never 
coalesced and stars never formed. We have light, but light 
in selective places. 

It has been pointed out that the abrupt cutoff of the 
luminosity function at the bright end and the flat slope 
at the faint end compared with the Press-Schechter mass 
function anticipated by the hierarchical clustering scenario 
suggests that there is a relative deficiency of light at high 
and low mass extremes compared with the situation at in- 
termediate masses (Yang et al. 2003, van den Bosch et al. 
2003). Semi-analytic models that follow the transforma- 
tion of gas into stars lead to similar expectations (Blanton 
et al. 1999, Somerville et al. 2001, Ostriker et al. 2003). 
These same results are found directly in observations. A 
full discussion of the observational situation is provided by 
Tully (2004). A condensed version is presented here. 



2 The High End of the Mass Spectrum 

Two lines of evidence strongly indicate that the ratio of 
blue light to mass decreases by about a factor 6 in progress- 
ing from halos of mass 1O 12 M to halos of mass 

2.1 Virial Analysis of Groups 

Collapsed halos are the sites of groups and clusters of 
galaxies (where in some circumstances the 'group' may be 
a single object). The relative distribution and motions of 
the galaxies in a group or cluster provide an indication of 
the mass of the ensemble. Many groups are relatively young 
(crossing times are a substantial fraction of the age of the 
universe) so are not expected to be relaxed. As a conse- 
quence, mass estimations based on the virial theorem have 
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Figure 1. Correlation between mass/light ratio and group 
crossing time. Left panel: all groups with 5 or more members in 
the volume limited T87 sample. Right panel: Subset of groups 
with at least 6 members with Mb < —17 distinguished ac- 
cording to whether a majority are E-SO-Sa (filled squares) or 
Sab-Sc-Irr (open circles). 



uncertainties of as much as a factor 2. Nonetheless, bear 
with this analysis based on the virial assumption because 
the variances in mass to light with environment that will 
be demonstrated are much larger than a factor 2. 

A volume limited sample of groups in the Local Super- 
cluster, extending to 25h^ 5 Mpc (/175 = Ho/75), was as- 
sembled by Tully (1987: T87), with the individual galaxies 
in the T87 groups identified in Table II of Tully (1988). The 
base galaxy catalog was considered complete to M B within 
this volume. The group catalog consists of 179 groups of 
two or more galaxies, 50 with five or more. 

The left panel of Figure 1 reveals a striking correlation 
between the virial mass My to blue light Mb ratio and 
the dynamical timescale of groups as measured by crossing 
time tx- There is an increase of an order of magnitude in 
My I Lb as tx decreases from ~ 0.5Ho to < 0.1 Ho. The 
right panel of the figure shows a subset of the sample, 
groups with at least 6 non-dwarf members, distinguished 
by whether the majority of the members are early or late 
type. There is a clear separation: the spiral rich groups have 
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Figure 2. Mass vs. light for groups over the full range of the ob- 
served mass spectrum. The large open circles and filled squares 
represent the same data seen in the right panel of Fig. 1. The 
crosses represent the remainder of the groups with 5 or more 
members found in the left panel of Fig. 1. The small open cir- 
cles represent groups with only 2-4 known members. The 3 filled 
inverted triangles represent E/SO/Sa groups more distant than 
the volume limited sample. The 5 filled triangles with error bars 
represent groups at the low end of the mass spectrum that are 
discussed separately. The My /Lg =94 dotted line is the mean 
value for the sample reported by T87. The solid line is a fit con- 
strained by the data with a high mass slope of 7 = 0.7 and an 
exponential fall off at masses less than M' = 6 X 1Q 11 Mq. 

long dynamical times and low My /Lb values compared 
with the groups dominated by E/SO galaxy types. 

Mass and luminosity are separated in Fig. 2. It is seen 
that the E/SO/Sa dominated groups tend to be the most 
massive. Two concurrent trends are found: more massive 
groups tend to have fractionally less blue light, and at a 
given mass E/SO/Sa groups have less blue light than their 
counterparts dominated by later morphological types. 

The properties of groups from the volume limited sur- 
vey region with only 2-4 known members are included in 
Fig. 2 in order to explore conditions in a lower mass range. 
In these cases the errors in the mass estimates are very 
large, which accounts for the large scatter. The 5 data 
points with error bars in this figure are discussed later. 



2.2 Numerical Action Methods 

Compatible results have emerged from the reconstruction 
of galaxy orbits with Numerical Action Methods (Peebles 
1989, Shaya et al. 1995, Tully and Shaya 1998). These re- 
constructions extend dynamical studies to larger scales. 
One is considering the dynamical influence of groups on 
each other and potentially taking into account mass that 
extends beyond the groups. For this discussion, no attempt 
is made to model the complex orbits within groups. Groups 
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Figure 3. Virgo infall constraints from two Numerical Action 
models. Data points indicate the velocities and separations from 
the center of the Virgo Cluster of individual galaxies. Galaxies 
within the 6° caustic of the cluster are indicated by open cir- 
cles. Galaxies outside the caustic but associated with the infall 
region are indicated by filled circles. The vertical brackets are 
located at angles from the center of the cluster that intersect 
infalling groups, so at lincs-of-sight that have received attention 
in the Numerical Action models. The amplitudes of the brack- 
ets illustrate the range of infall velocities anticipated by the two 
models under consideration. Left panel: M/Lb = 200 assigned 
to all units. Right panel: M/Lb = 1000 assigned to the Virgo 
Cluster and other E/S0 knots and M/Lb = 150 otherwise. 



are considered as ensemble halos and it is the orbits of the 
centers of mass of these ensembles that are reconstructed. 

In fact it is not the specific orbits that are recon- 
structed that we take seriously but rather it is the mass 
requirements that are found to give good fits. The recon- 
struction is constrained by the need to match the angular 
positions and radial velocities of objects in a catalog. If the 
cosmology is set with a specified age then the main vari- 
able is the mass assignments. With a given set of mass as- 
signments, orbits are found which leaves each entity today 
with 6 specified values in phase space, 3 of which rigorously 
match the observed boundary conditions. One of the other 
3 phase space values - the distance - is subject to test. Our 
model constraint is a \ 2 comparison of model and observed 
distances. The orbits that are found are not unique given 
the limitations of the catalogs and the poor constraints on 
the 3 free phase space components, but only models within 
a modest range of mass assignments give good \ 2 fits to 
the distance constraints. Basically, increasing mass forces 
objects with a given velocity to be farther removed from 
the main mass centers. 

In the earliest trials the minimalist assumption was 
made that M/Lb values are constant in all environ- 
ments. The subsequent best fits were found for M/Lb ~ 
2OOM /Lq . However this modest M/Lb assignment fails 
utterly to account for the large infall velocities toward the 
Virgo Cluster that are observed. Rather, a value M/Lb ~ 
IOOOMq / Lq is required for the cluster. Figure 3 attempts 
to summarize the observational constraints. The ampli- 
tude of the envelope of observed infall velocities is closely 
governed by the mass of the cluster. A very high mass of 
1.3 x 10 Mq is required to explain the velocity breadth of 
the infall region. Upon assigning the corresponding large 
mass to light ratio to the cluster (and other E/S0 dom- 
inated groups) it was found that the mass to light re- 
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quiremcnt for the remainder of the objects was reduced 
to M/ Lb ~ 15OM /L . 

The mass requirement found by the Numerical Action 
Method for the Virgo Cluster is 60% higher than found 
with the virial analysis, but applies to the cluster on a 
larger scale. The overall average M/Lb requirement of the 
action modelling is also higher than found via the virial 
theorem for most groups but again refers to larger scales. 
Both the virial and action studies suggest that M/Lb val- 
ues increase by a factor 5-7 in proceeding from the environ- 
ments of spiral galaxies to the Virgo Cluster environment. 



3 The Low End of the Mass Spectrum 

The faint end of the luminosity function of galaxies in 
low density environments has a substantially flatter slope 
than the theoretical mass function of hierarchical cluster- 
ing (Klypin et al. 1999, Moore et al. 1999). During a study 
of possible variations of the luminosity function with en- 
vironment (Tully et al. 2002) the idea arose that groups 
with low mass with very little light might be common. It 
was recalled that during the construction of the T87 group 
catalog a number of entities that were called 'associations' 
had group-like qualities except in one respect. These enti- 
ties had similar, though somewhat smaller, dimensions to 
loose groups and similar, actually smaller, velocity disper- 
sions. If considered to be bound, then the inferred masses 
are low - in the range 10 11_12 Mq. If one ignores the prop- 
erties of the galaxies and regards them only as test parti- 
cles probing a potential then one would conclude that these 
T87 associations are just the low mass end of a continuum 
of groups conditions. 

The outstanding difference with familiar groups with 
higher mass is that these low mass entities have very lit- 
tle light. The candidate members of these groups are all 
dwarf galaxies. Though the inferred masses are low, the 
luminosities are extremely low and Mj Lb values are high 
- so high that at the time of T87 there was a reluctance 
to conclude that the entities are bound (though there was 
a suspicion that this was the case!) 

Today, with greater boldness, it is argued that these 
groups of dwarfs are bound. It is to be noted that candidate 
groups of dwarfs are not rare. Within the restricted do- 
main of 5 Mpc, Tully et al. (2002) report 4 good examples. 
This number is comparable to the number of established 
groups within the same distance, and the census of dwarfs 
is surely incomplete at low Galactic latitudes. As discussed 
by Tully (2004), there is now a dramatic improvement in 
knowledge regarding the distances of nearby galaxies. Hub- 
ble Space Telescope observations resolves the brighter stars 
and permits a distance determination from the magnitude 
of the tip of the red giant branch (Karachentsev et al. 
2003). Good distances now exist for almost all the dwarf 
group candidates and it is found that the spatial correla- 
tions are equally as good in distance as they were found 
to be on the plane of the sky and in velocity. One of the 
groups of dwarfs turned out to be considerably closer than 
suspected based on velocities, with all the observed dis- 
tances concordant at the nearer position. 

At this point the evidence very strongly suggests that 
the groups of dwarfs are bound. If so, mass estimates sug- 



gest M/Lb ~ IOOOMq/Lq for these groups. The struc- 
tures are almost surely not virialized because crossing 
times are a substantial fraction of the Hubble time. Conse- 
quently, the derivation of masses are subject to large errors. 
The uncertainties can easily be a factor 2 and conceivably 
as bad as a factor 3. However, the M/Lb values that are 
found are an order of magnitude larger than the values 
attributed to loose groups of spirals. There is little doubt 
that if the groups of dwarfs are bound then they are reser- 
voirs of mostly dark matter. 

The mass and light values of the groups of dwarfs are 
shown in Fig. 2 as the filled triangles with error bars. In 
fact, five points are plotted in this manner and one stands 
apart in terms of luminosity. This higher luminosity entity 
is the near side component of what is commonly referred 
to as the Sculptor Group, but which is probably two sep- 
arate bound units. The near component includes NGC 55 
and NGC 300 plus smaller objects. The dynamical prop- 
erties (dimensions and velocity dispersion) suggest a very 
low mass for this proposed group, like that of the other 4 
groups of dwarfs. However, the two NGC galaxies are inter- 
mediate luminosity galaxies, not dwarfs, so the group is not 
deficient in light. It is to be concluded that not all group 
halos in the mass regime 1O 11_12 M lack light. However in 
the ensemble there seems to be a cutoff in light proceeding 
to low mass halos. 



4 The Dark Side 

The solid line in Fig. 2 is a fit of the equation 

L B = CM"'e- Mi/M (1) 

This equation relates luminosity and mass with 3 con- 
straints: a logarithmic slope at the high mass end, 7, a 
low mass exponential cutoff set by M', and a normaliza- 
tion, C. The fit minimizes a \ 2 expression with uncertain- 
ties in mass. Perhaps a more dramatic presentation of the 
same information is seen in Figure 4, a plot of M/Lb vs. 
mass (appreciating that the groups with 2-4 members rep- 
resented by the small open circles have very large uncer- 
tainties in mass). It is seen that halos in the mass range 
1O 12_13 M manifest the most light. Groups with progres- 
sively more mass up to 1O 15 M have progressively smaller 
light fractions, with this trend particularly enhanced in the 
groups dominated by early type galaxies. At the other end 
of the mass spectrum, below 1O 12 M there is an abrupt 
cut-off, with M/Lb values becoming very large and per- 
haps going off to infinity. One can suppose that there are 
low mass halos completely devoid of stars. 

The census of galaxies within 25h^ Mpc of T87 that 
lead to the group catalog can be used to construct a halo 
mass function. This is seen in the top panel of Figure 5. 
The mass function is poorly constrained at the high mass 
end by the small number statistics in such a restricted vol- 
ume. At the low mass end, the error bars only account for 
the y/N bin counting uncertainties. There are much larger 
horizontal errors associated with the uncertain masses of 
small groups. In the lowest mass bins the counts are dom- 
inated by individual galaxies outside of groups, in which 
cases it was assumed that M/Lb = 100. 
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Figure 4. M/Lb vs. mass over the full range of density 
regimes. The data are the same as in Fig. 2. The solid curve and 
dotted line are transpositions of the curve and line in Fig. 2. 



Admitting that there are considerable uncertainties 
in this mass function because of the restricted sample and 
low mass cutoff, perhaps the greater interest is in the lower 
panel of Fig. 5. Here one see the fraction of the total mass 
found in the separate mass bins. The filled part of the his- 
togram is contributed by the E/SO/Sa groups, while the 
open part represents the contribution by the groups domi- 
nated by spiral and irregular galaxies. One entity, the Virgo 
Cluster, contains ~ 40% of the mass in the survey volume 
(although it contains only 15% of the light!). The groups 
of early types together contain 60% of the mass. Consid- 
ering groups of all type, 90% of the mass is in groups with 
M > 3 x 1O 12 M . 

The groups of dwarf galaxies have only been discov- 
ered within a restricted 5 Mpc volume. Within this space, 
though they appear to be numerous, they only contain 
~ 10% of the mass associated with luminous galaxies. Yet 
the luminous groups within 5 Mpc lie on the tail of the his- 
togram in Fig. 5 just above 10 12 A/q. It is concluded that 
these groups of dwarfs, though they are dominated by dark 
matter, they make only a minor contribution to the total 
mass budget of collapsed dark matter halos. From a cos- 
mological perspective, it is the high mass to light ratios in 
massive clusters with short dynamical timescales that have 
greater implications. The force is with this dark side. 
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Figure 5. Top panel: Observed mass function in the volume 
within 25/i^g 1 Mpc and | 6 |> 30. The numbers of objects in 
each bin are given along the top of the lower panel. In brackets 
are the numbers of individual galaxies; ie, galaxies outside of 
groups. Lower panel: Fraction of the total mass associated with 
the mass function in half-dec bins. The mass associated with the 
groups of predominantly early types is indicated by the filled 
histograms. The open histograms correspond to the mass found 
in the groups with predominantly late types. 
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